Mott Insulators, No-Double-Occupancy, and Non-Abelian Superconductivity 
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SU (4) dynamical symmetry is shown to imply a no-double-occupancy constraint on the minimal 
symmetry description of antiferromagnetism and d-wave superconductivity. This implies a maximum 
doping fraction of \ for cuprates and provides a microscopic critique of the projected SO(5) model. 
We propose that SU(4) superconductors are representative of a class of compounds that we term 
non-abelian superconductors. We further suggest that non-abelian superconductors may exist having 
5(7(4) symmetry and therefore cuprate-like dynamics, but without d-wave hybridization. 
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Suppression of double occupancy on sites in the cop- 
per oxide planes is critical in explaining why cuprate sys- 
tems are antiferromagnetic Mott insulators at half filling 
and become superconductors through hole doping. The 
symmetric Zhang 50(5) model Q] predicts no charge 
gap at half filling. To recover Mott insulator phases 
at half-filling in the Zhang model it is normal to im- 
pose a no-double-occupancy rule by Gutzwiller projec- 
tion. This breaks 50(5) symmetry, but lattice calcula- 
tions and schematic arguments suggest that many 50(5) 
features might survive in such a projected 50(5) model 

We have proposed a unified description of high temper- 
ature superconductivity and antiferromagnetism based 
on a U(A) D SU(4) dynamical symmetry that has an- 
alytical solutions in three symmetry limits [f| 0. The 
50(4) limit of the 5J7(4) model corresponds to an anti- 
ferromagnetic phase, the SU (2) limit to a d-wave super- 
conducting phase, and the 50(5) limit to a critical sym- 
metry interpolating between the antiferromagnetic and 
superconducting phases. 

Although the methodology of the 5C/(4) model differs 
substantially from that of the Zhang model, its 50(5) 
limit represents the Zhang 50(5) algebra subject to con- 
straints implied by embedding 50(5) in the larger al- 
gebra SU(4). In this Letter we address the physical 
understanding of why SU (4) should play a crucial role 
in high temperature superconductivity, how no-double- 
occupancy and Mott insulator properties lie at the basis 
of this understanding, and provide a microscopic under- 
standing of the projected 50(5) model. 

The U (4) D SU (4) model has 16 symmetry generators: 
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Qij = E^^U^-fcj fti = (?«)* (!) 
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where cl i creates an electron of momentum k and spin 



projection i,j = 1 or 2 = f or I, Q = (jr, 7T, it) is an AF 
ordering vector, f2/2 is the electron-pair degeneracy, and 
g{k) is the d-wave form factor 

g(k) = (cos k x — cos k y ) » sgn(cos k x — cos k y ). (2) 

By forming new linear combinations, Q may be replaced 
by operators having more direct physical meaning: 
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ft = (^ t ) t , D^=p\ 2 , D= Pl2 , M=±(5 n +5 22 ) 

where Q+ is associated with charge density waves, 5 is 
the spin operator, Q is the staggered magnetization, D* 
(D) is the creation (annihilation) operator of spin-singlet 
d-wave pairs, 7P (7?) are associated with spin-triplet pairs, 
and 2M = n — f2 is related to the number (charge) oper- 
ator n. The representation space of the 5(7(4) model is 
built by the coherent D and 7r pairs: 

\n x n y n z n s ) = (4)^(4)^(^(0^0). (4) 

The operator Q + commutes with all generators and will 
be ignored in this discussion. Thus the most general ef- 
fective Hamiltonian in the symmetry-dictated truncated 
space is a linear combination of all scalar products con- 
structed from the remaining 15 5/7(4) generators. As- 
suming only one and two-body interactions, 
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where e, v, Go, Gi, x, and g define effective microscopic 
strengths of single-particle and interaction terms. 
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FIG. 1: (Color online) A schematic hole pair. Fuzzy balls 
are sites where electron holes form a pair: one hole at r, the 
other with equal probability (1/4) at the four neighboring 
sites (r = r ± a and rib). Balls connected by bright lines 
are sites where the presence of a hole would imply double 
occupancy. 



The operator set ® is closed under t/(4) D t/(l) x 
St/ (4) symmetry (hereafter termed SU (4)) only if the 
approximation in (J2J holds. The physics implied in this 
approximation becomes more transparent if we transform 
JTJ to coordinate space (using the exact form of Eq. 
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Hj = ^(-) r Cf jCr,* (6) 
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= c i,i c r,j — 
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where cj i {c r ,i) creates (annihilates) an electron of spin 

i located at r and c\ i (c F ,i) creates (annihilates) an elec- 
tron of spin i at the four neighboring sites, r±a and rib, 
with equal probabilities (a and b are lattice constants in 
x and y directions, respectively), 
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The factor (-) r in Eq. JSJ is (-)"*+"« and (n x ,n y ) are 
the coordinates of a lattice site on the copper oxide plane, 
b, which is positive for even sites (n x i 
n y = even) and negative for the odd sites [n x + n y = 
odd). This factor originates from the assumption e l ®' T w 
(— ) r and implies Mott insulator properties: the electrons 
are localized at lattice sites with small overlap between 
orbitals of electrons on neighboring lattice sites. 

From the coordinate representation JJJJ we see that 
spin-singlet and spin-triplet pairs are formed by holes on 
adjacent sites. Fig. 1 illustrates the spatial structure of a 



hole pair: if one hole is at r, the other hole occupies the 
four adjacent sites (ria and rib) with equal probability. 
The summation over r in the pair creation (annihilation) 
operators indicates that such pairs are highly coherent. 
It also can be seen that 
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with (n^) and (S^) the total electron number 
and spin operators at even (odd) sites, respectively: 
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Thus Q+ and Q represent the differences in total charge 
and spin between even and odd sites, respectively. 

However, @ does not close under commutation unless 



(12) 



(that is, c\ i (cf,i) is a basis for particles occupying sites 
adjacent to r). This separates lattice sites into categories 
A and B: if r = even are A sites (with operators cj i 

( c r.i)), r = °dd are the B sites (operators cj i (c F .i)), or 
vice versa. Then (11211 permits © to be written as 
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with Qij = i where the i sign is + (— ) if A 

is chosen to be even (odd) sites. (Whether A sites are 
taken to be even or odd is a labeling choice and does not 
influence the physics.) Then by explicit commutation the 
operators l|13f) close an SU (4) algebra. But by Eq. Q , 



\c r ',i >4,j} 
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jEs(*)<W+t (14) 



-1 for t = i2a, i2b, 

-1 for t = +a±b,-a±b 



(15) 



and H12(l is generally not satisfied unless the second term 
on the right side of (|14|l can be ignored. This term van- 
ishes if a constraint is imposed that whenever there is an 



3 



electron pair c^-c^ at r (see Fig. 1), no pair is permit- 
ted at r' = r + t, leaving nothing to be annihilated by 
Cf>i- This is a no-double-occupancy constraint because 
without it there is a finite amplitude for double site oc- 
cupancy. For instance, if one pair is at r' = r + 2a and a 
second pair at r, the probability is for two electrons to 
be located at r + a (see Fig. 1). We conclude that closure 
of the SU (4) algebra is a direct consequence of no double 
occupancy in the copper oxide conducting plane. 

Additional insight follows from observing that the va- 
lidity of (|12[) actually follows from the more general re- 
quirement that no pairs overlap, a consistency condition 
ensuring that the pair space and the pairing correlations 
be well defined, is sufficient to satisfy |Q . The no-pair- 
overlap constraint implies naturally that if a pair is cen- 
tered at r, no pair may be located at r' = r + 1 with t 
given m (US), and thus Eq. (JUJ) holds. 

For an A^-dimensional basis the minimum closed alge- 
bra is SO(2N) if all bilinear particle-hole and pair op- 
erators are taken as generators. The simplest basis for 
cuprates may be regarded as 4-dimensional since elec- 
trons can exist only in four basic states, on ^4-sites or 
.B-sites, with spin up or down. Thus, absent further con- 
straints, the minimum Lie algebra for the set of gener- 
ators that can describe high T c superconductivity and 
antiferromagnetism simultaneously in a cuprate system 
is 50(8) and not SU{A). The 28 generators of 50(8) are 
the 16 operators in (|13fl plus the 12 operators 
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(± depends on the even-odd choice for A sites; see H13H . 

Equation i|16|) contains two kinds of new (spin-singlet) 
pairs created by p\ 2 and q\ 2 , which may be termed 5 and 
5* pairs, respectively. In both the two electrons (holes) 
occupy the same site, with equal probability to appear 
anywhere in the lattice coherently. The 5* pairs differ 
from 5 pairs in their phases. The operators Sij are the 
hopping operators with and without spin flip, and Qij is 
the staggering of the hopping. These operators change 
D and tt pairs into 5 and 5*, or vice versa. 

The 50(8) algebra reduces to 5Z7(4) if the 5 and 5* 
pairs may be neglected, which occurs if we assume on- 
site Coulomb repulsion pushing the 5 and 5* pairs to 
sufficiently high energy. Thus, restriction to no double 
occupancy effectively allows the operators in Eq. (|16fl to 
be ignored and reduces 50(8) to the subalgebra 5f/(4). 
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FIG. 2: (Color online) Schematic spatial distribution for 
maximal hole pair occupation. Fuzzy balls are lattice sites 
where the electron holes form pairs. Balls connected by bright 
lines indicate sites where the presence of a hole would lead to 
double occupancy. 



Therefore, the minimal Lie algebra that can describe 
antiferromagnetism and d-wave superconductivity in a 
cuprate system is in general 50(8), but under the con- 
straint of no double occupancy the symmetry effectively 
reduces to 5f/(4). The assumption of an 5t/(4) symme- 
try in a cuprate system automatically implies the impo- 
sition of a no- double- occupancy constraint on the general 
50(8) symmetry in the copper-oxygen planes. 

It is likely that the hopping operator 5y in Eq. (|To|) is 
the source of the most important SU (4) symmetry break- 
ing terms. It breaks SU(4) but is a generator of 50(8), 
so this perturbation may be taken into account by an 
extension from SU(A) to SO (8). However, we may ex- 
pect the no-double-occupancy rule and thus the 5C/(4) 
symmetry to be a good initial approximation. 

The implicit 5t/(4) occupancy constraint dictates an 
upper limit for the doping fraction in SU (4)-conserving 
states. Fig. 2 illustrates the spatial distribution when 
the number of hole pairs is maximal. By counting, the 
maximum number of holes is Q = O e /4, where Q e is the 
total number of lattice sites. Thus the largest doping 
fraction preserving SU(4) symmetry is Pf = Q/Q e = j. 
The empirical maximum doping fraction (0.23 ~ 0.27) 
for cuprate superconductivity may then be interpreted 
as a direct consequence of 5t/(4) symmetry. 

The preceding discussion implies that: (1) The physi- 
cal origin of cuprate SU (4) symmetry is proximity of an- 
tiferromagnetism and d-wave pairing, coupled with sup- 
pression of double occupancy by on-site Coulomb repul- 
sion. (2) Two important facts in cuprates, that normal 
states are Mott insulators and that superconductivity ex- 
ists only in a narrow doping range (P < 0.27), are direct 
consequences of an SU (4) dynamical symmetry. 

Superconductivity in cuprates is a specific example of 
what we shall term non-abelian superconductivity, which 
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differs from conventional superconductivity in the rich- 
ness of pair structure for condensed states and in the ap- 
pearance of competing sources of long-range order. The 
key issues for 5C/(4) non-abelian superconductivity are 
that coherent pairs are formed by holes on adjacent sites 
so that both singlet and triplet states contribute, and 
that alternative long-range order (antifcrromagnetism) 
enters on an equal footing with superconductivity. 

In contrast to BCS superconductivity, which is de- 
scribed by a single dynamical symmetry chain having 
only abelian subgroups (5/7(2) D U(l)), the minimal 
symmetry consistent with cuprate data is 5/7(4), which 
has a much richer structure (three dynamical symmetries 
having non-abelian subgroups and differing fundamen- 
tally in their properties). We propose that the differences 
in observational characteristics for these two types of su- 
perconductivity originate in this difference in dynamical 
symmetry structure and in non-abelian superconductiv- 
ity resulting from electron-electron interactions instead 
of electron-phonon interactions. 

The primal role of 5/7(4) symmetry in non-abelian 
cuprate superconductivity suggests that any pairing 
structure leading to the 5/7(4) algebra entails dynamics 
similar to that of cuprates. Therefore, cZ-wave symmetry 
of the pairs need not be critical to non-abelian supercon- 
ductivity in general and 5/7 (4) superconductivity in par- 
ticular. Pairs with any internal symmetry (extended s- 
wave, p-wave, mixed symmetry, . . . ) could exhibit 5/7 (4) 
superconductivity if the no-double-occupancy constraint 
is valid and correlations can form adjacent-site pairs. 
Generally, c\ i may be defined as 

4,i = £i. 9 (t)i 2 = i (17) 

t t 

where t is a few finite lattice displacements of r and git") 
is the form factor. Different forms of g(t) reflect different 
internal symmetries of the pairs, but they all satisfy the 
condition (T^jl under no double occupancy and thus pre- 
serve the 5/7(4) algebra and the general Hamiltonians 
implied by its dynamical symmetry chains. The struc- 
ture J7J) of the d-wave pairs is only a special case of l(T7j) 
with t = ±a, ±b and ff(±a) = \ and g(±b) = — |. 

In summary, we have shown that 5/7(4) is the min- 
imal symmetry accommodating superconductivity, anti- 
ferromagnetism, and a no-double-occupancy constraint 
in cuprate systems, and that 5/7(4) symmetry implies 
a maximum doping fraction of \ in the cuprates, by 
symmetry alone. Because the Zhang 5*0(5) algebra is 
a subalgebra of 5/7(4), these results indicate that closure 
of the Zhang algebra also implies no double occupancy. 
Why then does the SO (5) model require Gutzwiller pro- 
jection? The work presented here suggests that the pro- 
jection requirement arises from assumptions inconsistent 
with the underlying 5*0(5) algebra in the Zhang effective 
Lagrangian formulation, where five of the 577(4) genera- 



tors (D^,D,Q) are treated as order parameters forming 
a superspin vector. Thus one cannot apply algebraic con- 
straints to them through the commutators. This is most 
easily seen if not only the 50(5) generators but also the 
elements of the Zhang superspin vector are treated as 
operators rather than order parameters, thereby enlarg- 
ing the algebra to 5C/(4). Within the 5/7(4) framework, 
there is no need for projection. As demonstrated in Refs. 
0,0], within the parent 5/7(4) group antiferromagnetism 
and superconductivity are described by different dynam- 
ical symmetries (50(4) and 5/7(2), respectively), and 
50(5) is a critical dynamical symmetry that interpo- 
lates between SO (A) and £17(2). In the 5*7(4) model 
it is the 50(4) symmetry, not the 50(5) symmetry, that 
naturally describes undoped states, and the spectrum for 
unbroken 50(4) dynamical symmetry is intrinsically an- 
tiferromagnetic with gapped charge excitations. 

Our results imply some important consequences of at- 
tributing cuprate superconductor behavior to an 5C/(4) 
algebra that follow directly from symmetry, independent 
of details: (1) Normal states are antiferromagnetic Mott 
insulators. (2) Hole doping of normal compounds leads 
first to 50(5) fluctuations in both antiferromagnetic and 
superconducting order (implying phases that may exhibit 
spin glass or stripe character), and then to 5/7(4) non- 
abelian superconductivity. (3) 5/7(4) superconductivity 
is strongly suppressed for doping fractions exceeding ^. 
(4) Symmetry breaking resulting from violation of the 
no-double-occupancy constraint may be described by a 
parent 50(8) algebra where terms that break SU (4) may 
still respect 50(8) symmetry. (5) 5/7(4) symmetry, not 
d-wave pairing per se, is the ultimate cause of cuprate be- 
havior, implying that systems could exist having non-<i 
pairing but cuprate-like dynamics. The first three con- 
sequences are postdictions in strong accord with exist- 
ing data. The fourth is a prediction that may be tested 
through detailed applications of the SU (4) model to data. 
The final prediction may be tested by searching experi- 
mentally for compounds having pairing structures other 
than d x 2_ y 2 that satisfy the 5/7(4) algebra. 

We thank Aditya Joshi for discussions and for careful 
checking of the 5/7(4) commutation relations. 
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